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Collisions are a major modification process over the history of the Kuiper Belt. Recent work illuminates the 
complex array of possible outcomes of individual collisions onto porous, volatile bodies. The cumulative 
effects of such collisions on the surface features, composition, and internal structure of Kuiper Belt Objects 
are not yet known. In this chapter, we present the current state of knowledge of the physics of cratering and 
disruptive collisions in KBO analog materials. We summarize the evidence for a rich collisional history in 
the Kuiper Belt and present the range possible physical modifications on individual objects. The question of 
how well present day bodies represent primordial planetesimals can be addressed through future studies of the 
coupled physical and collisional evolution of Kuiper Belt Objects. 



1. INTRODUCTION 

The Kuiper Belt contains some of the least modified ma- 
terial in the solar system. Some Kuiper Belt objects (KBOs) 
may be similar to the planetesimals that accreted into the 
larger bodies in the outer solar system. However, KBOs 
have suffered various modification processes over the life- 
time of the solar system, including damage from cosmic 
rays and ultraviolet radiation, sputtering and erosion by gas 
and dust in the interstellar medium, and mutual collisions 



(e.g. Stern 2003 1. Robust interpretations of the surfaces and 



internal structures of KBOs require improved insight into 
the relative weight of each of these processes. 

The present understanding of the importance of colli- 
sions on the physical evolution of KBOs is limited by the 
state of knowledge in three fundamental areas: (1) the dy- 
namical history of the different populations within the trans- 
Neptunian region (Morbidelli et al. this volume); (2) the 
physical properties of KBOs {Brown, Stansberry et al. this 
volume); and (3) how the physical properties of KBOs (ex- 
pected to be icy and porous) affect the outcome of colli- 
sions (this chapter). The dynamical history of a population 
defines the evolution of mean impact parameters (velocity, 
angle, mass ratio of the projectile and the target) within 
and between KBO populations. The impact parameters and 
the material properties of the colliding bodies determine the 
outcome of an individual impact event. Finally, the cumu- 
lative effects of collisions are determined by the coupled 
physical and dynamical evolution of KBOs. 

Variable progress has been made in these three areas. 
Over the past decade, great improvements in observations 
and models have illuminated the rich dynamical history of 
the Kuiper Belt. At present, there is a sparse, but growing, 



body of data on the physical properties of KBOs (e.g., size, 
density, composition, and internal structure). Although a 
significant body of work has been devoted to collisions be- 
tween icy, porous bodies, our understanding of the govern- 
ing physics is still incomplete. The collisional evolution of 
KBOs is a particularly interesting and challenging problem 
because of the range of possible outcomes that depend on 
the changing dynamical structure of the Kuiper Belt. 

In this chapter, we present a summary of the work to 
date that can be applied to the physical effects of collisions 
in the Kuiper Belt. We begin with observational evidence 
for significant past and present-day collisions in the Kuiper 
Belt (Q. We then present a range of possible outcomes 
from collisions between KBOs (Q and discuss the prin- 
cipal discriminating factors (Sj4j). Based on the expected 
physical properties of KBOs, we summarize the results of 
laboratory and numerical experiments that have been con- 
ducted to determine how material properties, such as com- 
position, porosity and impact conditions, including velocity 
and mass ratio, affect collision outcomes (Sj5). Finally, we 
discuss several open questions and future research direc- 
tions for studying collisions in the Kuiper Belt. (f|6]). 

2. EVIDENCE FOR A RICH COLLISIONAL HIS- 
TORY IN THE KUIPER BELT 

In this section, we summarize four observations that sup- 
port a significant collisional history within the Kuiper Belt. 
First, we discuss observations of interplanetary dust parti- 



cles (IDPs) by the Pioneer and Voyager spacecraft ([2.1 



Analyses of the orbits of IDPs conclude that the Kuiper Belt 
must be one of the dust source regions. Second, the size dis- 
tribution of Kuiper Belt Objects has at least one break from 
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a simple power law around diameters of 10's km, which 
is consistent with models of collisional equilibrium among 
the smaller bodies ({2.2 \. Third, the discovery of a possi- 



ble dynamical family of objects in the Kuiper Belt implies 
conditions that produced at least one near-catastrophic col- 



lision of one of the largest KBOs ({23 I. Finally, models of 
the accretion of the largest KBOs demonstrate that the mass 
in the ancient Kuiper Belt must have been much larger than 
observed today. The total mass loss of >90%, and perhaps 
as much as 99.9%, was driven by a combination of dynam- 



ical perturbations and collisional grinding ({ 2.4 1. 

In addition to the observable features discussed below, 
collisions within a small body population will also affect 
rotation rates, surface colors, and the formation of bina- 
ries. The rotation rates of bodies in collisional equilibrium 
will reflect the angular momentum transfer from typical im- 
pact conditions (see e.g., Love and Ahrens\\991 Paolicchi 



et al. |2002 for asteroid rotations). The formation of binary 
KBOs is still a matter of debate. Some binaries seem to 
have formed via collisions, while others have too much an- 
gular momentum for a collision origin (Margot 2002). The 
observed color diversity in the Kuiper Belt is also contro- 
versial and not correlated directly with collision energy (see 
Doressoundiram et al. this volume). However, the range of 
outcomes from collisions depend on material properties as 
well as the impact parameters. The growing data on rotation 
rates, colors, and binaries will provide in the future addi- 
tional constraints on the collisional evolution in the Kuiper 
Belt. 

2.1. Interplanetary Dust Particles 

Dust and small meteoroids were detected in the outer 
solar system by the Pioneer 10 & 11 and Voyager 1 & 2 
missions (Humes 1980, Gur nett et q/.||199"7) l. Pioneer 10 
and 1 1 measured the concentration and orbital properties of 
dust from 1 to 18 AU. The dust impacts detected by Pio- 
neer 1 1 between 4 and 5 AU were determined to have either 
high inclination or eccentricity or both. In other words, the 
IDPs were either not on circular orbits and/or not on near 
planar orbits. Hence, the observed increase in particle flux 
at Jupiter could not be explained by gravitational focusing, 
which is inefficient for highly inclined and eccentric orbits, 



and Humes ( 1980) > suggested that the dust had a cometary 
origin. 

In a reanalysis of the Pioneer data, Landgraf et al. { 2002 ) 
found that the dust flux was relatively constant at distances 
exterior to Jupiter's orbit. To produce a constant dust flux 
from drag forces, the dust must originate from a source be- 
yond the detection locations by the spacecraft. \Landgrqf 



et al. ( 2002 1 modeled the radial dust contribution using 
three source reservoirs, dust from evaporating Oort Cloud 
and Jupiter family comets and dust from collisions between 
KBOs, and argue that the amount of dust observed by Pio- 
neer 10 and 11 can only be explained by a combination of 
all three reservoirs. They find that comets can account for 
the material detected inside Saturn's orbit but an additional 



reservoir is necessary for the dust observed outside Saturn's 
orbit. 

Although Voyager 1 and 2 did not carry specialized de- 
tectors for dust, \Gurnett et al. ( \ 991) found that the plasma 
wave instruments could detect impacts from small particles 
with masses > 10~ u g (two to three orders of magnitude 
below the mass detection limit by Pioneer 10 and 1 1). From 
data collected between 6 and 60 AU, \Gurnett et q/.| ( [T997] l 
found a severe drop off in dust detection events after 5 1 AU 
and 33 AU for Voyager 1 and 2, respectively. As a result, 
the authors conclude that the source of the dust cannot be 
interstellar. Furthermore, the small latitudinal gradient de- 
creases the likelihood that the source objects are planets, 
moons, or asteroids (if the dust did originate from such ob- 
jects, one would expect a strong latitudinal gradient since 
the planets, moons, and asteroids are effectively all in the 
same plane). The Voyager IDP observations are consistent 



with a dust source from the Kuiper Belt ( Gurnett et al. 1997 
\Jewitt and Lmm|2000| i. 

In summary, the radial distribution and orbital properties 
of outer solar system IDPs cannot be explained by source 
material solely from Jupiter family comets and Oort cloud 
comets and indicate the need for an additional active source 
of dust in the outer solar system. The IDP observations 
are well matched by models of dust produced during the 



collisional evolution of the Kuiper Belt (e.g.,\Jewitt andLuu 



2000 and i 2.4 1. Dust derived from mutual collisions in the 



present day Kuiper Belt is analogous to the zodiacal dust 
from the asteroid belt (Mtiller et al. 2005 ) and observations 
of rings of dust around other main sequence stars (Moro- 
Martin et al. and Liou et al. this volume). Because the 
removal time of dust is much shorter than the age of the 
solar system (Kenyon et al. this volume), the dust must be 
replenished by collisions occurring throughout the history 
of the solar system. 

2.2. Size Distribution of KBOs 

Formation models indicate that KBOs accreted within a 
thin disk with low relative velocities and inclinations (Mor- 
bidelli et al. this volume). However, the present velocity 
dispersion (~ 1 km s" 1 ) and the inclination distribution 
(about 20 degree half width) of KBOs are both much higher 



than expected during the coagulation stage \Trujillo et al. 



2001 ). The large relative velocities and the large inclination 
distribution of the KBOs point to significant dynamical in- 
teractions with Neptune, which resulted in a rich collisional 
history \Davis and Farinella\\991) . 

If the bodies in the Kuiper Belt were fully collisionally 
evolved and collision outcomes were independent of size, 
the differential size distribution (dN ~ r~ q dr, where N is 
number of bodies in the size bin of radius r) would be de- 
scribed by a self-similar collisional cascade and fit by a sin- 
gle power law index of q = 3.5 ( Dohnany t"|1969 Williams 
and Wetherill 1994). If the population is only partially col- 
lisionally evolved and/or the disruption criteria is depen- 
dent on scale, the size distribution will deviate from a single 
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power law. For example, the size distribution in the aster- 
oid belt deviates from a simple power law in part because 
of strength effects {O'Brien and Greenberg 2005) and re- 
cent collisions, such as dynamical family forming events 
^dell'Oro etal.\200\) . 

Recent observations indicate that the size distribution 
of KBOs has a break at diameters of 10's km, with fewer 
smaller bodies than expected from extrapolation from bod- 
ies of 100's km diameter (Bernstein et al. 2004, Chen et a~L\ 
2006 Roques et al. 2006). The slope of the differential size 
distribution of large KBOs (< 25 magnitude, >100 km di- 
ameter) is well established, with a slope in the range of 4 to 



4.8 ^Trujillo et alp OOl Petit et al. 2006, Petit et al. this 
volume). Bernstein et al. ( 2004[ ) also suggested that the 
classical KBOs have a different size distribution from the 
other dynamical populations (for KBO population classifi- 
cations, see Gladman et al. this volume). 

Over the past decade, several groups have made sig- 
nificant progress in modeling the collisional evolution of 
the Kuiper Belt [ Davis and Farinella\\991\ Stern and Col 



well\l997~\\Kenyon and Bromley\2004\ Pan and Sari\2005 
Kenyon and Luu 1999] Kenyon et al. this volume). Their 
work provides a theoretical basis for a break in the KBO 
size distribution around 10's km. \Davis and Farinella\ 
\991\ first demonstrated that few of the largest bodies in 



the Kuiper Belt experience catastrophic disruption events 
in which 50% of the mass is permanently removed. In other 
words, most of the largest KBOs are primordial; they have 
persisted since the end of the coagulation stage, although 
some may have suffered shattering collisions. 

Collision evolution models indicate that the break in the 
size distribution corresponds to the upper size limit of the 



collisionally evolved population (Davis and Farinella 1997 
Kenyon and Br omley 2 004 Pan and San||2005y Over 
time, collisions preferentially disrupt smaller objects be- 
cause of their lower critical disruption energies and their 
higher number densities (and hence higher collision prob- 
abilities) compared to larger bodies. When the disruption 
criteria is size-dependent, the collisionally evolved size dis- 
tribution deviates from q — 3.5 (see O'Brien and Gree~n~\ 
berg\2 QQ3 ). For example, Pan and Sari ( 2005 \ utilize a dis- 
ruption criteria proportional to the gravitational binding en- 
ergy of the body, and the equilibrium power law has q = 3. 
Numerical evolution simulatio ns by \Kenyon and Bromley\ 
( 2004 ) and analytical work by \Pan and Sari\ ( |2005| ) are in 
good agreement with the observation of the number of 10's 
km size bodies by \Bernstein et al. (2004). Note that the lo- 



cation of the break in size between the collisionally evolved 
and primordial populations increases with time and depends 
on the dynamical evolution of the system. 

The size distribution of KBOs is likely to have a sec- 
ond break in slope at significantly smaller sizes. The sec- 
ond break corresponds to the transition between the colli- 
sionally evolved strength-dominated bodies and collision- 
ally evolved gravity-dominated bodies. As the strength of 
KBOs is essentially unknown at present, the location of the 
strength to gravity transition based on catastrophic disrup- 



tion models, e.g., 10's to 100's m (S 5.2 1, is highly uncertain. 

In summary, the observed size distribution of KBOs de- 
parts from a single power law, which is consistent with the 
existence of both a collisionally evolved population and a 
primordial population ( Davis and FarineUa\\991\ Kenyon 



and Bromleyl2004\ \Pan and S ari 2005; Stern and Colwell 



1997 ). In this scenario, the primordial populations should 



have experienced primarily surface modification processes 
through impact cratering events (S|3] Figs. [T] and |2^b). The 
population of largest bodies probably has a subpopula- 
tion with differentiated internal structures (Mer k and Pri-\ 
alnik 2006 ) and a subpopulation with rubble pile structures 
(Davis and Farinella 1997). Both the collisionally evolved 
strength and gravity-dominated populations would have ex- 
perienced the full range of collisional outcomes, including 
catastrophic disruption and changes in internal structure and 
composition (fj3| Fig. |2j;). 

2.3. Kuiper Belt Family (2003 EL61) 

At present, a few tens of dynamical families have been 
identified in the asteroid belt (Bendjoya and Zappala 2002 ). 
These objects are grouped together in proper element space 
and have similar spectral features where detailed observa- 
tions are available. The orbits of the family members can be 
integrated back in time to a common starting point, suggest- 
ing formation via a catastrophic collision. Although colli- 



sional evolution models of the Kuiper Belt (?2.2i indicate 



that very few KBOs larger than about 100 km have experi- 
enced catastrophic disruption events, Brown et al. (2007b ) 
have observed what seems to be a collisional family in the 
Kuiper Belt (see also Brown this volume). 

2003 EL61 has two known satellites and five proposed 
family members. All of these objects have similar proper 



elements, colors, and a deep H2O spectral feature (Brown 
\et a/.|2007b| ). Although the detection of H2O ice on the sur- 
faces of KBOs is not unique to these objects, the significant 
depth of the spectral feature is characteristic of the proposed 
family, suggesting either more recent or more abundant ex- 
posure of surface ice compared to other KBOs. 2003 EL61 
has a double-peaked rotational light curve with a period of 
3.9 hours (Rabinow itz et q/.|2006| ). Assuming that the light 
curve is solely due to the equilibrium shape of a rotating, 



homogeneous, fluid ellipsoid, Rabinowitz et al. (2006]) and 
\Lacerda and J ewitt (2007]) derive the size (~1500 km di- 
ameter) and density (~ 2.6 g cm -3 ). However, the derived 
size and density are highly uncertain as 2003 EL61 is likely 
to possess nonzero shear strength (Holsapple 2007 ). 

These combined observations suggest that 2003 EL61 
suffered a significant, but sub-catastrophic impact event 
( Brown et q/.|2007"b ). If the modeled bulk density of 2.6 g 



cm" 3 " is correct and the pre-collision density of 2003 EL61 
was comparable to other large KBOs (~2 g cm -3 , Brown 
this volume), then about 20% of the original mass was lost. 
In this model, the dispersed material was preferentially H2O 
ice, presumably derived from an ice-rich mantle, produc- 
ing the shared water spectral feature of the proposed family 
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members. 

Further investigation of the proposed 2003 EL6 1 family 
and search for other dynamical families would provide use- 
ful constraints on the collisional history of the Kuiper Belt. 

2.4. Total Mass of the Kuiper Belt 

The total mass in the modern Kuiper Belt is depleted 
from a smooth surface density extrapolation from the gi- 
ant planet region of the solar system. Based on the ob- 
served size distribution of bodies between 30 to 50 AU, 
the total mass is only about 0.01 Earth masses (less than 
5 Pluto masses) ^Bernstein et a/.||2004| l. However, |fffem| 



( 1996} and \Stern and Colwell ( 1997| l demonstrate that the 
Kuiper Belt must have been more massive in the past for 
the largest KBOs (100 to 1000 km) to accrete via mutual 
collisions. At least 90% of the mass in the Kuiper Belt was 
lost through collisions and ejections induced by the stirring 
and migrating of Neptune ^Stern and Colweil\\991 Hahn 
\and Malhotra\\999) . 

3. POSSIBLE COLLISION OUTCOMES IN THE 
KUIPER BELT 

The observations summarized in the previous section in- 
dicate that collisions are an important factor in the evolution 
of the Kuiper Belt. In studying the physical effects of col- 
lisions between KBOs, we are guided by the mature stud- 
ies of collisions in the asteroid belt { Asphaug et a/.|[2002 



\Holsapple et a/.|2002 i. However, the possible outcomes of 
collisions between KBOs are more diverse compared to as- 
teroids because of the dynamical state of the system and the 
range of physical properties of individual KBOs. 

The important observations that inform the range of pos- 
sible collision outcomes are as follows. The dynamical 
state of the Kuiper Belt has changed dramatically with time 
(Morbidelli et al, this volume); hence the mutual colli- 
sion velocities between KBOs also varied with time. In 
the modern Kuiper Belt, the mutual collision velocities are 
around 1 km s _1 { Trujillo et a/.||2001| l for classical KBOs 
and slightly higher for the other populations (Gladman et 
al. this volume). For these impact velocities, most bodies 
smaller than 100's km size have experienced a catastrophic 
disruption event, while most of the larger bodies have sur- 
vived. All bodies should have suffered the production of a 
significant density of surface impact craters. 

In addition to the dynamical impact conditions, the phys- 
ical properties of KBOs are important. While the present 
observations are limited, the range of bulk densities of 
KBOs is <1 to ~ 2.6 g cm -3 (Stansberry et al. this vol- 
ume) and the largest KBOs have a range of surface volatile 
compositions (H2O, CH4, etc) in addition to a refractory 
(rock and organic) component (Barucci et al. this volume). 
From these observations and studies of short-period comets, 
believed to be fragments from KBOs, a typical KBO has 
a significant (but unknown) fraction of volatiles and high 
porosity. 

Given the possible range of material properties and im- 
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Fig. 1 . — Top: Schematic of an impact crater on a target made 
of a mixture of volatile and refractory material. The energy of 
the impact produces melting and devolatization at the base of the 
crater and in the ejecta. Bottom: Schematic of an impact crater on 
a target made mostly of volatile material. The surface of the target 
is covered with a "crust" of darker refractory material. The impact 
is large enough to excavate fresh volatiles from depth creating a 
bright ejecta blanket. Collapse of the crater wall creates a darker 
region at the bottom of the crater. 
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Fig. 2. — Schematic showing a possible evolutionary track for a 
KBO. a) A cross section through a primordial planetesimal with 
low bulk density and high microporosity. b) After many small im- 
pacts, the bulk density has increased and volatile composition de- 
creased at the surface, c) After a large sub-catastrophic collision, 
the body is shattered and the surface is covered with ejecta. d) A 
catastrophic impact event disrupts the body, creating a rubble pile 
with high macroporosity and heterogeneous internal composition. 
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pact conditions, we outline the potential array of dramati- 
cally different outcomes from collisions between KBOs: 

1. Surface Impact Features: As a result of the rich 
collisional history of the Kuiper Belt, impact craters 
are expected to be common on the surfaces of KBOs. 
However, the morphologies and size distributions are 
sensitive to the surface and internal structure of the 
body. Some highly porous bodies survive the for- 
mation of multiple, large craters comparable to the 
radius of the object, as on the low-density asteroid 
Mathilde \Veverka et aZ.||1997| l. Alternatively, only 
small craters may be observed on rubble piles formed 
by catastrophic disruption, such as asteroid Itokawa 
( Fujiwara et a/.|| 2006). 



2. Surface Composition and Color: Impact craters 
and catastrophic disruption events may darken the 
surface by removing volatiles via heating from the 
energy of the impact or brighten the surface of a body 
by excavating fresh ices (Fig. [TJ. 

3. Density of Surface Materials: Laboratory impact 
craters in highly porous and compressible materials 
compact the impact site, creating density heterogene- 
ity on the surface \Housen et a/.||1999] l. Over time, 
cumulative small impacts on a microporous surface 
may increase the bulk density and decrease the bulk 
porosity (Fig. |2^b). In contrast, modeling results in- 
dicate that nearly all ejecta from a crater in a macro- 
porous body may reach escape velocities, leaving the 
bulk density unchanged (Asphaug et al. 1998 ). 

4. Internal Structure and Composition: A sub- 
catastrophic impact may shatter a body (and create 
a large crater) but leave the original internal material 
relationships intact (Fig. |2j;), while a catastrophic im- 
pact both shatters and disperses a body such that the 
gravitationally reaccumulated remnants are rubble 
piles with high macroporosity and mixed composi- 
tion (Fig. |3Ji). 

In the next section, we describe the factors that control 
the outcome from individual collision events. 

4. FACTORS THAT CONTROL OUTCOMES OF 
COLLISIONS 

Recent advances in the understanding of the physics of 
collisional processes between icy, porous bodies provide 
new fuel to the study of the role of collisions in the Kuiper 
Belt. The outcome of collision events are governed by the 
impact conditions (velocity, angle, and mass of each body) 
and the physical properties of the colliding bodies (strength, 
composition, and internal structure). Both the impact con- 
ditions and physical properties affect the efficiency with 
which the energy of the impact is coupled to the target. 
In this section, we summarize four overarching factors that 
control the outcome of an impact event between KBOs. In 



the following section (f|5]), we will describe laboratory and 
numerical experiments on KBO analog materials that inves- 
tigate these controlling factors. 

First, the composition and internal structure of the bod- 
ies determines the critical velocity required to enter the 
strong shock regime, where the deformation and coupling 
of energy and momentum can be described through the 
Rankine-Hugoniot conservations equations (S4.1 i. Slower 



impact events, where plastic deformation dominates, re- 
quire more detailed knowledge of the physical properties 
(particularly strength) of the bodies compared to the strong 
shock regime. Collisions between KBOs are likely to span 
the range of plastic and shock deformation. 

Second, the final outcome is a balance between the 
forces of strength and gravity 4.2 1. Scaling laws have 



been developed for cratering and catastrophic disruption in 
each regime, but a large transition region exists. Because of 
the low gravity and expected low strength of KBOs, many 
collisions may fall in the transition region. 

Third, the internal structure and composition of the col- 
liding bodies may significantly affect collision outcomes 



(E 4.3 1. Some of the impact energy will be partitioned into 
phase changes when highly volatile materials are present. 
High levels of porosity also alter the energy coupling by 
acting as a shock absorber and localizing shock deforma- 
tion. The momentum coupling with high porosity changes 
the excavation flow in the cratering regime and the disper- 
sal of fragments in the disruption regime compared to col- 
lisions between solid bodies. 

Fourth, collision outcomes are sensitive to the mass ratio 
of the projectile and target. At the same kinetic energy, a 
larger projectile is more efficient at removing mass than a 



smaller projectile (S 4.4 1 



4.1. Shock Deformation 

In most high energy impact events, the deformation is 
driven by a shock wave. The energy from the shock con- 
trols the physical deformation from the collision, such as 
fragmentation, pore collapse, heating, and phase changes. 
The shock also determines the deposition of momentum 
that leads to crater excavation or dispersal of fragments fol- 
lowing a catastrophic disruption event. The amount of de- 
formation can be estimated by considering the volume of 
material shocked to a given peak pressure. 

A strong shock wave is produced in a hypervelocity im- 
pact event, where the impact velocity exceeds the bulk 
sound speed (q,) in both the target and projectile. How- 
ever, mutual collision velocities between KBOs are likely 
to span the range from subsonic to supersonic (hyperveloc- 
ity) collisions (see below). When collisions are comparable 
to the sound speed, plastic deformation dominates, rather 
than strong shock deformation. Under subsonic conditions, 
collisions are simply elastic and governed by the coefficient 
of restitution. 

In this section, we present a summary of the shock 
physics that determines the outcome in the hyperveloc- 
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ity regime. The peak shock pressure is deduced from the 
conservation equations and material equation of state, de- 
scribing the pressure-volume-temperature (P — V — T) 
states. A shock wave satisfies the Rankine-Hugoniot (R-H) 



and the equations for nonporous and porous H2O ice are 
given in 1 



Stewart and Ahrens (2004 2005 ). 



mass, momentum, and energy conservation equations {Rice 
[19581 : 



Ui - u 
Ei — Eq 



Us ( 1 

Us 
Vq 
1 
2 



Vq 

(Ui - Uq) , 

(Pi - P Q ) (Vq - Vi) 



(1) 

(2) 
(3) 



In the above formulae u is particle velocity, Us is shock 
velocity, V is specific volume (= l/p, where p is density), 
P is pressure, and E is specific internal energy. The initial 
unshocked state is subscripted and the final shocked state 
is subscripted j. 

The shock Hugoniot is the curve that describes the locus 
of possible P — V shock states for a given initial P — V — T 
state. For a given impact scenario, the shock pressure is cal- 
culated using Eq. |2j the impact velocity, and the equations 
of state of the target and projectile. Many materials may be 
described using a simple linear Us — u shock equation of 
state of the form $Ruoff\1967) , 



Us 



c + su, 



(4) 



where c and s are material constants (for their relationship 
to finite strain theory, see Jeanloz\1989 1. The linear shock 
equation of state is simply a representation of the P — V 
shock Hugoniot translated into Us — u space using the R-H 
equations. 

In the planar impact approximation (also called the 
impedance match solution, see derivation in Melosh 1989), 
the particle velocities induced by the shock wave reduces 
the projectile's velocity and mobilizes the target such that 
continuity at the projectile-target interface is achieved and 



(5) 



where v is the impact velocity and subscripts t and p refer to 
the target and projectile, respectively. The shock pressure is 
derived by solving for u t using the equality of Eq. [2] in the 
target and projectile and substituting Eqs.[4]and[5] 

Po.t( c t + s t u t )u t = po,p(c p + s p (v - u t ))(v - u t ). (6) 

The quadratic function for ut is readily solved. In the case 
of identical shock equations of state in the target and projec- 
tile, the particle velocity is equal to v/2, and the peak pres- 
sure is given by p (c+sv/2)(v/2). Because of strength and 
phase changes, the Us — u shock equations of state for natu- 
ral materials are usually fit with multiple linear segments in 
u, corresponding to different pressure ranges on the shock 
Hugoniot. The shock equations of state for many rocks and 
minerals are compiled in \Ahrens and Johnson ( 1995a|b| l, 



As the shock wave propagates into the target, the peak 
pressure, derived from the planar impact approximation, 
decays from rarefaction waves on the free surfaces. The 
size of the region at peak pressure (known as the isobaric 
core) and the decay exponent depend on the impact velocity 
and material properties (e.g., equation of state and poros- 
ity) ^Ahrens and Q'Keefe\\l9S7\ \Pierazzo et~aT\\991\ . In 
general, the pressure decay is steeper for high velocities be- 
cause of energy partitioning into phase changes. The occur- 
rence of impact-induced phase changes can be estimated by 
considering the critical shock pressures required for melting 
and vaporization. When the shock pressure is above a crit- 
ical value, the material is melted/vaporized after passage of 
the shock wave and return to ambient pressure conditions. 

The present mean mutual collision velocity between 
classical KBOs (~1 km s _1 ) is lower than the bulk sound 
speed of full density silicates and ices. Nonporous H 2 
ice has a q, of 3.0 km s _1 at 100 K (Petrenko and Whit- 



\worth\\999\\Stewart and Ahrens\2QQ5). Silicate ro cks have 
larger q,, typically around 5 km s _1 (Poirier 2000 1. Sounds 



speeds of laboratory preparations of nonporous ice-silicate 
mixtures, with up to 30% by weight sand, are similar to pure 
H 2 ice (Lange and Ahrens 1983). Pure porous H 2 ice, 
on the other hand, can have much lower sound speeds, from 
0.1 to 1.0 km s _1 for bulk densities of 0.2 to 0.5 g cm -3 
{Mel lor\\l975\ \Furnish "a nd Remo 1997). Silica aerogels 
with densities of about 0.2 g cm have sounds speeds of 
about 200 m s _1 [ Gross et a/.|l988 i, and 35% porous sand 
has a sound speed of 130 m s _1 . If KBOs are volatile rich 
and porous, then mean present-day collisions may be super- 
sonic. During the collisional evolution of the Kuiper Belt, 
collisions span the subsonic to supersonic regimes. 

Understanding the controlling physics in the subsonic 
regime, where plastic deformation dominates, will require 
focused studies on analogs for the range of mechanical 
structures in the Kuiper Belt (Sj5j. In the strong shock 
regime, crater scaling relationships and catastrophic disrup- 
tion theory are applicable, as described in the next section. 

4.2. Strength and Gravity 

The final outcome of a collision, e.g., crater size or dis- 
persed mass, depends on the balance between strength and 
gravitational forces. In the case of impact cratering, the 
relationship between the size and velocity distribution of 
the impacting population and the observed crater popula- 
tion can provide insight into the collisional history of a sys- 



tem (as has been done for the terrestrial planets [Strom et al. 



|2005| ) and asteroids ^O'Brien and G reenberg 2005 1). Back- 
ing out the impactor properties requires the application of 
the appropriate crater scaling relationships, which depend 
on both the impact conditions and material properties. In 
the case of catastrophic disruption, knowledge of the prop- 
erties of the populations of disrupted and primordial bodies 
provide strong constraints on the collisional evolution of the 
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Kuiper Belt. Here, we present the crater size and catas- 
trophic disruption scaling laws in the strength and gravity 
regimes. 

4.2.7. Crater Scaling Theory 
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Fig. 3. — Cratering efficiencies (ratio of ejected and displaced 
mass to projectile mass, nv) as a function of the ratio of gravita- 
tional to inertial forces (inverse Froude number, 7r2) for different 
target materials at different velocities. Lines denote fitted crater- 
ing efficiencies in liquid water, dry sand (35% porosity), and weak 
rock using Eq.[l0] Experimental data (see text): nonporous ice A : 
50% porous ice •; various crushable non-icy mixtures with porosi- 
ties of about 40% □, 70% O, and 96% V- Colors denote impact 
velocity. 

Development and validation of the appropriate scaling 
relationships is crucial for the next generation of collisional 
evolution models of KBOs that include consideration of 
physical deformation effects. Because of their low gravity 
and likely low strength, the outcome of collisions between 
KBOs is near the transition between the strength and gravity 
regimes. In this section, we discuss the strength to gravity 
transition and summarize the equations and material param- 
eters for crater scaling for comparison to laboratory craters 



in ice and porous targets in ( 5. 1 



In a cratering event, the shock-driven excavation flow 
produces a roughly hemispherical cavity, called the tran- 
sient crater. Assuming that material strength can be repre- 
sented by a single parameter, Y, the transition size between 
the strength and gravity-controlled cratering regimes is pro- 
portional to Y/pg. Y is the dominant strength measure that 
controls crater size (e.g., shear strength); p is the bulk den- 



sity of the target; and g is the gravity of the target §Melosh 
1977 Melosh and McKinnon\1978 \. As the impact energy 
increases, the outcome of collisions will transition from a 
cratering regime to a total body disruption regime. The 
criteria for catastrophic disruption, Q* D , is defined as the 



specific energy (kinetic energy of the projectile divided by 
the mass of the target) required to disrupt and gravitation- 
ally disperse half the mass of the target Melosh and Ryan 



1997). Note that, unlike the cratering regime, disruption is 
governed by the bulk tensile strength of the body, which is 
typically an order of magnitude lower than the compressive 



strength of brittle solids (see §4.2.2| &|44 1. 

The theory for crater size scaling based on impact pa- 
rameters and material properties is summarized bylTfoZsap"- 



ple (1993). A common approach utilizes 7r-scaling, with 
empirical constants derived from impact and explosion cra- 
tering experiments under Earth's gravity and high gravity. 
Predicting the final crater volume and shape requires two 
steps: (i) calculating the volume of the transient crater cav- 
ity using the 7r-scaling laws and (ii) calculating the amount 
of collapse of the transient crater to the final crater volume 
and shape. The first step is better understood than the sec- 
ond. 

In 7r-scaling, the cratering efficiency, 7Ty, is defined as 
the ratio of the mass of material ejected and displaced from 
the transient crater cavity to the mass of the projectile: 



7Ty 



pV 



Mc 



(7) 



where V is the volume of displaced and ejected material, 
m p is the mass of the projectile, and M c = pV. For 
strength-dominated craters, the cratering efficiency depends 
on the ratio of a measure of the shear strength of the target, 
Y, to the initial dynamic pressure, given by 



7Ty- 



Y 



(8) 



where, v± — v sin 9, v is the impact velocity, and 6 is the 
impact angle from the horizontal. In the gravity-dominated 
regime, the cratering efficiency depends on the ratio of the 
lithostatic pressure at a characteristic depth of one projectile 
radii, r p , to the normal component of the initial dynamic 
pressure (the inverse Froude number): 



T2 



2 ' 

± 



v 



(9) 



where g is the gravitational acceleration. 

Impact experiments demonstrate that the transition from 
strength to gravity-dominated cratering spans about two 
decades in 7r 2 . Following \Holsapple (1993]) and \Holsapple 
and Housen ( 2004|), the cratering efficiency can be defined 



by an empirical, smoothed function of the form 

k(tt 2 



71V 



P/a 
7T - 
Y 



(10) 



where the exponents are related to a single coupling expo- 
nent, p, by a — 3p/(2 + p) and (3 — 3p/2. The coupling 
exponent p is bounded by two cratering regimes: momen- 
tum scaling (where p = 1/3) and energy scaling (p = 2/3) 
^Holsapple and Sch midt 1987; Holsapple 1987). Note that 
Eq. [10] assumes that the target and projectile have the same 
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density. The transition from strength to gravity dominated 
regimes occurs when Y ~ pgr p . 

In Fig. [3] cratering efficiencies are presented for liquid 
water (K=0.9_8, fi = 0.55, Y = MPa), dry sand (K=0.132, 
/i = 0.41, Y = MPa, 35% porosity), and weak rocks 
(K=0.095, p = 0.55, Y = 3 MPa) (values from |tfo/^qp^fe 
i/oMsew|2004| ). Dry sand is a non-crushable porous ma- 
terial and weak rock is a reasonable analog for nonporous 
H 2 ice. Cratering efficiencies in crushable, porous materi- 
als, from hypervelocity experiments in vacuum, lie a factor 
of few below the dry sand line [Schultz et al. 2005). The 
transition from strength regime (lower values of 772, when 
Try > 7T2) to gravity regime (higher values of 7^) corre- 
sponds to the transition from a horizontal line in a iry — 7r 2 
plot, when the cratering efficiency is independent of tt2, to 
a power law with slope —a. The cratering efficiency in the 
strength regime increases with increasing impact velocities, 
as indicated by the curves for impacts into weak rock tar- 
gets at 0.5, 2.0, and 7.0 km s _1 . Note that the cratering 
efficiency in dry sand is less than for weak rock in the grav- 
ity regime because of energy dissipation in the porous sand. 
Data from impact cratering experiments conducted in vac- 



uum under Earth's gravity into nonporous ice (A, \Cintala 



et al. 



2005 



1985; Lange and Ahrens\\92>l\\Burchell and Johnson 



and 50% porous ice (•, Koschny et aZ.|2001 ; Burchell 



et al. 2002} fall in the strength regime. Cratering exeri- 
ments at 1.86 km s^ 1 in porous mixtures of sand and per- 
lite bonded with fly ash and water under one atmosphere of 
pressure and varying gravity are nearly independent of 7T2, 
indicating strength-dominated behavior with a plausible in- 
tersection with the gravity-dominated regime (the dry sand 
Housen and Ho lsapple 2003| l. 



line) (§5.1.3 



Cratering events in the Kuiper Belt by a nominal 0.5 m 
radius body at 1 km s^ 1 onto targets of 0. 1 to 1000 km radii 
correspond to 7r 2 values in the range from 10 -10 to 10 -6 . 
Therefore, for the range of impact velocities in the Kuiper 
Belt, Fig. [3] demonstrates that the presence of any strength 
is likely to control the final crater size for the majority of 
impact events. In the upper range of possible values of 7T2, 
gravity may dominate if cratering is less efficient in KBOs 
than in dry sand. 

After formation of the transient crater cavity by the ex- 
cavation flow, most craters undergo collapse to a final crater 
shape (Melos h and Ivanov\ [l999). For simple, strength- 
dominated craters, the final crater size is similar to the tran- 
sient cavity with some collapse of the crater walls. For 
cratering in soils and rocks, the rim radius o f the tran- 
sient crater is approximately R r = 1.73F 1 / 3 [Holsapple 



T993J. Complex, gravity-dominated craters undergo signif- 



icant collapse of the transient cavity, and the final crater rim 
radius scales with the transient crater rim radius and grav- 
ity by fi comp iex( cm) = 0.37i?,.(cm) 1 086 ( 5 / WO ' 086 
lHolsapple\i993\ . 

The final state of crater formation in porous, icy bod- 
ies is not well understood. In Sj5] we discuss some of the 
laboratory experiments that provide our best guesses at the 
appropriate crater size scaling laws for the Kuiper Belt. 



4.2.2. Catastrophic Disruption Theory 

Based on the models of the collisional evolution of the 
Kuiper Belt, it is probable that a large fraction of bodies 
have suffered both cratering events and disruptive colli- 
sions. The catastrophic disruption criteria, Q* D , is the ratio 
of the projectile's kinetic energy to the mass of the target 
required to disrupt and disperse half the mass of the target. 
The criteria has two components ( Davis et al.\1979\ , 



Q*S + Qb, 



(11) 



where Q* s is the strength of the body to shattering and Qb 
is the gravitational binding energy of the target. 

The catastrophic disruption criteria for a rocky body is 
shown by the thick solid line in Fig.p](§5.2.2|and \Benz and 



Asphaug 1999). The critical energy is averaged over all im- 
pact angles. A head-on collision is most efficient, requiring 
about an order of magnitude less energy compared to the an- 
gle average \Benz and Asphaug\\999 Leinhardt et al. 2000 



Leinhardt and Richardson 2002 ). In the strength regime, 



where Q* s dominates, the critical energy decreases with in- 
creasing target size because tensile strength, the controlling 



strength measure, is scale dependent {Housen and Holsap 
pie 1999) . The larger the body, the larger the number of 
pre-existing natural flaws and the lower the tensile strength. 
In the gravity regime, pressure from the self gravity of the 
object increases the strength, following the shattered rock 
curve for head-on impacts [Melosh and Ryan 1997[ ). In 
the gravity regime, the gravitational dispersal criteria dom- 
inates over shattering by orders of magnitude. Note that the 
standard disruption criteria curves assume that the size of 



the projectile is small compared to the target (see f 4.4 1. 

The manner in which volatile content and porosity affect 
the disruption criteria is not well understood. Here, we es- 
timate the effects of each using nonporous and porous H2O 
ice as an example. There has been little work on catas- 
trophic disruption of large objects in the gravity regime at 
impact speeds and compositions that are relevant to the the 
present day Kuiper Belt, thus this discussion is meant to 
provide general guidance not detailed values. 

In the strength regime, the Q* D intercepts for nonporous 
and porous ice (thin solid line and dashed line, respec- 
tively) are tied to results from laboratory disruption exper- 
iments (see §5.2. l| and Arak awa et a/.|2 002). These values 
are consistent with other experimental results ( Ryan et al. 
1999). The slope of Q* s for pure ice is assumed to be the 
same as for rock. For porous ice, on the other hand, the 
slope in the strength regime is particularly uncertain. The 
slope for porous materials may be much shallower than for 
nonporous materials because the size-dependent scaling of 
flaws may not apply \Housen and Hol sapple 1999). This 
uncertainty is depicted in Fig. |4]by several dashed lines of 
varying slope. Perhaps counterintuitively, in the strength 
regime, a porous material is harder to disrupt than a non- 
porous material due to localization of energy by compaction 
of pores and/or reflection of the shock wave off of free sur- 
faces (see §5.2. \\ . 
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Fig. 4. — Catastrophic disruption and dispersal energy (Qo) ver- 
sus target radius for rock (thick solid line), nonporous ice (thin 
solid line), 50% porous ice (dashed line). The negatively sloped 
portions of the curves are in the strength regime, the positive 
slopes are in the gravity regime. The criteria for rock is bas ed on 



angle averaged results from 3 km s collisions onto basalt i Benz 



and Asphaug 19991. The nonporous and porous ice intercepts are 
based on laboratory experiments (Arakawa et al. 2002 , QJ, = 40 
J kg -1 and Q* D = 143 J kg -1 for low and high porosity targets, 
respectively, of 5 cm radius). The extrapolation into the gravity 
regime is highly uncertain for porous materials. These results as- 
sume energy coupling by a small projectiles compared to the size 
of the target. 



In the gravity regime, the Q* D criteria for nonporous ice 
lies below the rock criteria by the ratio in mass (for this plot 
vs. target size). This is consistent with the numerical impact 
simulations into ice targets by \Benz and Asphaug ( 1999| l. 



Similarly, a porous target of the same size is easier to dis- 
rupt because of its lower total mass. Adjusting the gravity- 
dominated Q* D criteria by the ratio in total mass makes the 
unreliable assumption that the energy coupling from the 
collision is similar for each material. Because of the sig- 
nificant dissipative effects of porosity, porosity may have a 
large affect on energy coupling in the gravity regime (the 
Q* s term may be more significant). More work is needed to 
determine exactly how porosity affects the energy coupling 
for catastrophic disruption events. 

As with impact cratering events, it is clear that in order 
to predict the collision outcome from a disruption event, it 
is necessary to know something about the material prop- 
erties of the KBOs. In the case of catastrophic disruption 
global properties, rather than surface properties, are more 
important: for example, is the target porous, icy, rocky? In 



[5.2 we discuss the small amount of work on catastrophic 
disruption of KBO analog materials. 



4.3. Internal Structure and Composition 

KBOs are likely to possess a wide variety of internal 
structures, as depicted in Fig. [2] Dynamical excitation and 
increased collision frequencies from the migration of Nep- 
tune removed most of the mass from the original Kuiper 
Belt, leaving a mixture of collision fragments and unmodi- 
fied material ^Davis and Farinella 1997; Hahn and Malho-\ 
|fra|1999| Morbidelli et al. this volume). Comets may pro- 
vide clues to the present internal structure of KBOs; how- 
ever, comets are expected to be diverse themselves (Barucci 
et al. this volume, \Weissman et a/.|2 005 ). 

Porosity, either primordial or the result of collision 
events, is a major complicating factor in predicting the 
amount of shock deformation. Since KBOs are expected 
to contain a range of porosities, the outcome of individ- 
ual collisions could vary widely depending on the prop- 
erties of the colliding bodies. When the initial porosity is 
high, the shock impedance (the bulk density times the sound 
speed) is low, and the peak shock pressures produced for a 
given impact condition are lower compared to a solid target 
(Equations [T]& |2j. For a given shock pressure, however, 
the internal energy increase is larger in a porous material 
because of the greater change in volume during shock com- 
paction (Equation [3}. Hence, the temperature rise due to a 
shock is higher in porous materials, and impacts into porous 
ices may result in abundant melting or vaporization near the 
impact site {Stewart and Ahrens 2005, 2004). Porosity is 
an efficient dissipator of shock energy. As the shock prop- 
agates into the target, porosity increases the decay rate of 
the shock because of the increased energy partitioning into 
heat (e.g., Meyers 2001| l. Therefore, the shock-deformed 
volume in porous materials is smaller compared to a solid. 

The length scale of the porosity is also important. Small- 
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scale porosity compared to the shock thickness is described 
as microporosity . The thickness of the shock wave is pro- 
portional to the scale of the topography on the surface of the 
projectile. Large-scale porosity, e.g. a rubble pile of solid 
pieces, is described as macroporosity. In the latter case, the 
solid (e.g., monolithic) pieces may have high strength, and 
impact cratering events onto a monolithic piece would re- 
flect the high surface strength. For catastrophic disruption 
events, however, a rubble pile has low bulk tensile strength. 
In a rubble pile, the shock wave would reflect upon encoun- 
tering void space between solid boulders, and as a result, 
the energy from the shock would be deposited in a smaller 
volume compared to a shock wave propagating through a 
monolith of competent rock. On the other hand, a micro- 
porous body may have low surface compressive strength, 
but because of the efficient shock dissipation, a more ener- 
getic impact is required to catastrophic ally disrupt the body 
{Asphaug et a/.|[T998l. Hence, both macroporous and mi- 



croporous bodies may have high disruption strength. 

Compositional variation and surface layers also change 
the way energy is coupled into the target. The impact en- 
ergy will be partitioned into more compressible phases and 
a larger Q* D is required to disrupt a more compressible ma- 
terial {Benz and Asphau g\\999) . Because some of the en- 
ergy of the impact is partitioned into heating, each collision 
event will also result in net devolatilization. Finally, phase 
changes (melting, vaporization) of volatile materials will re- 
sult in steeper decay of the shock wave tyAhrens a nd Okeefe 
\T9Tl\\Pierazzo et al.\\99l\\Pierazzo and Melosh\2000) that 
tends to localize the shock deformation in a manner similar 
to the effects of porosity. 

4.4. Mass Ratio 

The mass ratio of the colliding bodies is also an impor- 
tant factor in determining the collision outcome. In sim- 
ulations of subsonic collisions, Leinha rdt and Richardson\ 
(2002 ) found that smaller projectiles were not as efficient at 
disrupting targets as larger projectiles (see Fig.Blfrom Lein 



hard t and Richardson\2062\ \Melosh and Ryan\1997ty . The 
mass ratio affects the volume over which the impact energy 
and momentum are deposited. When the projectile is much 
smaller than the target, the impact directly affects a small 
volume, about the size of the projectile. The rest of the 
target acts to dampen any material motion. When the pro- 
jectile mass is similar to the target mass, on the other hand, 
the projectile comes in direct contact with a significant vol- 
ume fraction of the target. As a result, the specific energy 
required to catastrophically disrupt a target decreases by or- 
ders of magnitude. 

The dependence on mass ratio has not been studied di- 
rectly for hypervelocity impacts, although the dependence 
on the coupling of energy and momentum should be sim- 
ilar to the subsonic case. First, the size of the peak pres- 
sure region (isobaric core) is proportional to the size of the 
projectile. Second, the decay of the peak shock pressure 
with distance depends on the impact velocity. The decay is 
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Fig. 5. — Mass of the largest reaccumulated body (M pr i m in 
units of total system mass M — M pro j + M ta rg) following a 
large collision event as a function of the kinetic energy of the pro- 
jectile. Results from iV-body impact simulations of km-sized rub- 
ble piles in the subsonic regime {Leinhardt and Richardson 2002). 
The crosses, filled hexagons, and open squares denote projectile 
to target mass ratios of 1:3, 1:6, and 1:9, respectively. In all cases, 
targets were identical with a radius of 1 km. All data points are av- 
eraged over several simulations at various impact angles (rms error 
bars). Note that, for the same kinetic energy, the largest projectiles 
produced the smallest largest post-collision remnant. These results 
indicate that increasing the projectile size increases the disruption 
efficiency. 
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steeper for higher impact velocities because more energy is 
partitioned into phase changes and deformation. Low im- 
pact velocities have a more shallow decay exponent (in the 
elastic limit). The particle velocities are proportional to the 
peak shock pressure; thus the shock pressure profile in the 
target will affect the dispersal of fragments and the catas- 
trophic disruption criteria. Q* D should decrease as the pro- 
jectile size increases for a fixed kinetic energy of the pro- 
jectile. 

Let us now consider a likely impact scenario in the 
Kuiper Belt. For an example 100 km volatile rich target 
in the Kuiper Belt, Fig. |i] predicts that about 10 5 J kg -1 
is necessary for catastrophic disruption. At 1 km s _1 , the 
projectile would have a radius two thirds that of the target. 
However, Fig. [4] assumes small (point source) projectiles. 
As the projectile to target mass ratio approaches unity, the 
amount of energy per target mass needed to disrupt the tar- 



surface layers of different strength materials can also have 



get drops (Holsapple 1993; 


Melosh and Ryan\\991\ Lein- 


hardt and Richardson\\2002 


i. Therefore, a 100 km target 



may indeed be catastrophically disrupted by a smaller pro- 
jectile than predicted in Fig. [4] and the larger objects in the 
Kuiper Belt may have suffered more catastrophic or near 
catastrophic impacts than inferred in previous studies. More 
work is needed to determine how Q* D behaves with mass ra- 
tio in hypervelocity collisions. 

5. STUDIES OF COLLISIONS IN ANALOGS TO 
KUIPER BELT OBJECTS 

We now turn to laboratory and numerical experiments 
on the major factors that affect the collision outcomes de- 
scribed above. Laboratory experiments and numerical mod- 
els of collisions between icy and/or porous bodies serve as 
the best analogs at present for impacts into KBOs. There is 
a large body of laboratory work on cratering impacts into 



volatile ices and mixtures as well as porous material (?5.1 
Although a coherent theory will require additional experi- 
ments, the laboratory results are a good guide to the possible 
outcomes of cratering collisions on KBOs. In comparison, 
the laboratory and numerical experiments on catastrophic 



disruption of KBO analogs are more limited (f 5.2 1. This is 
due in part to the inability to study catastrophic disruption 
in the gravity regime in the laboratory and the difficulties in 
modeling collisions between porous, volatile bodies. 

5.1. Cratering 

We begin with the results of several impact cratering 
studies into nonporous H2O ice and ice-silicate mixtures 
Then, the effects of porosity are introduced 
However, it is difficult to deconvolve the effects 



C'- 



5.1.1 



05.1.2 



of porosity and low material strength in laboratory studies. 
Possible outcomes include cratering events that result in 
compaction rather than the normal crater excavation flow 



(J 5.1.3 1. Relatively little work has been conducted on ices 
more volatile than H 2 0, which have been observed on the 
surfaces of the largest KBOs (Barucci et al. and Brown 



a significant affect on the crater morphology (t 5.1.5 1. Be- 
cause of the influence of an atmosphere on the final crater 



form, explosion cratering studies (Holsapple and Housen 



2004) are not included in this discussion. 

5.1.1. Cratering in nonporous H2O ice and ice-silicate 
mixtures 

As a low density and volatile material, solid H2O ice 
represents a very simple model for the bulk properties of 
KBOs. Depending on the evolution of KBOs, some sur- 
faces may be dominated by solid ice. Most laboratory 
impact experiments in ice are conducted in the strength 
regime. Generally, cratering efficiencies in solid ice are 
similar to a dry soil or weak rock (Fig. [3]) (Chapman and 



McKinnon\l986 I. Solid ice cratering experiments span im- 
pact velocities of 0.1 to 7.3 km s _1 using a wide range 



of projectile materials \Burchell and Johnson\\2005 Croft 


et al. 


1979| \Lange and Ahrens\\19S2\ 


1987 


\Shrine et al. 


2002 


Grey and Burchell 2003 


Kawakami et al. 1983 Kato 


et al. 


1995J Iijima et al\\\995 


Cintala et a/.||1985 1. For a 



given impact energy, the crater volume is more than an or- 
der of magnitude larger than craters formed in a typical hard 
silicate rock. In many of these experiments, the measured 
volumes of the craters include a component of spalled ma- 
terial (near-surface material ejected under tensile failure), 
forming a terraced crater morphology with a central pit. 
Hence, the reported volumes are larger than the transient 
crater volume and comparisons to 7r-scaling laws must be 
made with caution. Also, differences in ice target prepara- 
tion contribute to scatter between experiments. 

Impact cratering experiments in solid H2O ice have in- 
vestigated the effects of target temperature on the crater- 



ing efficiency \Lange and Ahrens\\\9%2\ [1987, Gre y and 
\Burch ell 2003). Low-temperature ice has a cratering effi- 
ciency between temperate ice and hard rock. Under fixed 
impact conditions, the crater depth and volume decreased 
by factors of 2 and 4, respectively, as the ice temperature 
decreased from 253 K to 100 K \Grey and~B urchell 2 003| >. 
It is well established that the yield strength of ice increases 
as the temperature decreases (e.g., Sammonds et al. | 1998 1; 
however, the magnitude of the effect is not well predicted 
\Grey and Burch ell 2003 ). 

Thermodynamic analyses of shock wave experiments on 
solid H2O ice at 100 K indicate that peak shock pressures 
of 1.6 and 4.1 GPa are required to produce incipient and 



complete shock-induced melting, respectively ( Stewart and 
Ahrens 2005). For pure ice on ice impacts, these pres 



sures are achieved at impact velocities of about 1 and 2 
km s _1 . If the bulk shock impedance of solid ice is sim- 
ilar to porous volatile-refractory mixtures, mass melting of 
solid ice within KBOs is only expected at the upper end of 
the range of collision velocities within the Kuiper Belt (0.5 
to 3 km s^ 1 \Dell'Oro et al.\200l\ . 



this volume) (i 5.1.4). Finally, experiments into targets with 



Nonporous mixtures of ice and silicates (e.g., ice- 
saturated sand) have also been studied in the strength 
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regime \Croft et al\\919 Koschny and Gnw||2001 1. The 
cratering efficiency decreases with increasing silicate con- 
tent. 

5.7.2. Cratering in porous material 

Several groups have conducted cratering experiments 



into porous materials (Schultz et al. 2005 ; Schultz and Gault 



[1985 ; Schmidt 19801 \Housen et al.\\999\\HblIsen and Hoi 



|^qp/?/e|[2~003| \Koschny et a/.pOOT I. The porous targets in 



elude porous ice, sand, Ottawa flint shot, pumice and vermi 
culite. Hypervelocity impact experiments, conducted under 
vacuum, into pumice powder with porosities between 35 
and 50% follow a single gravity-controlled crater scaling 
that is slightly less efficient than dry sand (Fig. ^ {Schultz, 
et al. 2005 1. However, for lower velocity impacts, in the 



strength regime or in the transition between strength and 
gravity-dominated cratering, the effects of porosity can be 



significant ({5.1.31. 

Results from cratering experiments into ~50% porous 
( Koschny et al. 2001 Burchell et al. 2002[ ) and solid ice tar- 
gets indicate that the displaced and ejected mass scales lin- 
early with impact energy. In other words, the crater volume 
is proportionally larger by the difference in target density. 
However, it is unlikely that this result can be extrapolated 
to events with much larger impact energies because of the 
considerable effects of vaporization on the final crater size 
(Holsapple and Housen 2007 ; Schultz et al. 2005 ). 

Because of the large increase in internal energy associ- 
ated with shock compaction of porous H2O ice, the criti- 
cal pressures required for shock-induced melting are lower 



compared to solid ice. From shock wave experiments, Stew 



art andAhrens (2004l find that shock pressures of only 0.3 



to 0.5 GPa initiate melting in 40-45% porous ice, and com- 
plete melting is reached by 2 GPa. These pressures corre- 
spond to impact velocities in the range of 1 to 2.5 km s _1 for 
collisions between porous ice bodies. Although porous ice 
has lower shock impedance than solid ice, the increase in 
internal energy from pore compaction results in similar crit- 
ical impact velocities for shock-induced melting. If KBOs 
have shock impedances greater than pure porous H2O ice, 
as expected if they are ice-refractory mixtures, then mutual 
collisions under the present dynamical environment will re- 
sult in abundant melting of H2O ice. 

Therefore, shock-induced melting in porous targets may 
produce crater cavities lined with quenched melt (rapidly 
cooled liquid that solidifies as a glass). Quenched melt 
lined craters have been observed in laboratory impact ex- 
periments into 50% porous H2O ice (nylon projectiles at 0.9 
to 3.8 km s" 1 ) ^Koschny et aZ.|2001 i and 5 to 60% porous 
soda lime glass (glass projectiles at 4.9 to 6.1 km s _1 {Love 
\et q/.|1993] l. In some cases, all the impact-generated melt 
was ejected from the crater. Hence, cratering events onto 
porous KBOs may produce solid ice ejecta fragments. 

The depth of penetration of the projectile plays a signif- 
icant role in the cratering efficiency in porous materials. In 
hypervelocity impacts, the impact angle determines the out- 



come. For example, the cratering efficiency in compress- 
ible porous perlite granules (p — 0.2 g cm -3 ) increases as 
the impact angle decreases from 90° to 30° ( Schultz et al. 
2005). Vertical impacts into porous materials penetrate 
deeply into the target, resembling a deeply buried explo- 
sion. Low angle impacts, in contrast, reach a more optimal 
shallow depth of burial to produce a larger crater. In the low 
velocity regime, an impedance mismatch between the target 
and projectile will also influence the depth of penetration. A 
dense projectile may experience little deformation and pen- 
etrate deeply, resulting in less efficient cratering compared 
to a projectile with density that matches the target. Interpre- 
tation of the cratering record on KBOs will need to include 
the role of impact angle and the depth of penetration in the 
final crater size. 

In the case where the projectile is more dense than the 
target (e.g., a solid rock meteoroid impacting a porous 
KBO), the impact conditions may be supersonic for the tar- 
get but subsonic for the projectile. In this case, the projectile 
is not significantly disrupted by the impact event. Labora- 
tory experiments show that the penetration depth increases 
as the density contrast between the projectile and target in- 
creases ( |Love et a7jl 993 ). Intact or melted nylon (p — 1.14 
g cm -3 ) and copper (p = 8.92 g cm -3 ) projectiles were re- 
covered after impacts at velocities up to 7 km s _1 into 50% 
porous H2O ice ( Koschny et al. |200T Burchell et q/.|2002] l. 
The experimental results suggest that dense meteoroids may 
embed themselves into the surfaces of KBOs and comets. 
In a pathalogical example, a population of compacted, de- 
volatilized projectiles might be found embedded in the sur- 
faces of very porous, volatile KBOs. 

Reliable numerical models of crater formation in porous 
materials have been hindered by the difficulty in modeling 
the shock compaction of porous materials (e.g., Herrmann 



1969 Johnson\\99\ I. Some general results can be drawn 
from the relatively few simulations to date: (i) a proxy 
model for porosity using layers of solid ice and void and 
the Autodyne code (Burchell and Johnson 2005 ), (ii) a P- 
alpha crush up model for sand using the CTH code $Housen\ 
and Holsapple 2000 ), and (iii) a new e-alpha compaction 



model using the iSALE code ]Wunnemann et a/.|2006| l. In 
summary, the transient crater diameters in porous materials 
are smaller but the crater is deeper than those in nonporous 
media. The lower bulk density of the porous target allows 
the projectile to penetrate more deeply. The shock wave is 
attenuated more quickly in porous material because energy 
is partitioned into crushing pores. These numerical experi- 
ments show that porous crushable objects are more resilient 
to large impact than nonporous objects because the damage 
from the impact is much more localized. With these more 
advanced models of porosity, future work can address the 
volume of material that experiences deformation (fragmen- 
tation, devolatilization) from impact events in the Kuiper 
Belt. 
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5.1.3. Compaction Cratering 

Observations of an unusual main belt asteroid, 253 
Mathilde, have incited several studies on the role of poros- 
ity on impact cratering. Imaged during a flyby of the NEAR 
spacecraft, Mathilde has a low bulk density (~ 1.3 g cm -3 ) 
and exhibits four large impact craters with diameters larger 
than the mean diameter of the asteroid ( Veve rka et al. | 1 997) . 
The large craters have no visible ejecta blankets or raised 
rims. As a result of their size, the craters are very close to 
each other and yet seem to show no evidence of interaction. 
The unique characteristics of Mathilde suggest that the in- 
ternal structure of this C-type asteroid is different from 
other classes of main belt asteroids in a fundamental way. 

\Housen et al. | (1999) and \Housen and H olsapple ( |2003[ ) 
conducted a series of cratering experiments into compress- 
ible, porous material in an attempt to explain the craters 
on asteroid 253 Mathilde. The authors suggest that high 
microporosity (40-60%) and high compressibility lead to a 
phenomena they termed compaction cratering. 

In their studies, the projectile and impact velocity was 
held constant, and the target porosity and gravity (using a 
centrifuge) were varied. In the high gravity environment, 
the craters had no raised rims and minimal ejecta outside the 
crater because most of the ejecta never escaped the crater 
cavity. A computed tomography scan of the crater showed a 
region of pore compaction approximately one crater radius 
below the crater. \Housen et a/.| ( |1999] ) and House n andHol-\ 
sapple ( 2003| > also impacted one of the used targets close to 
the original crater and confirmed that there was little inter- 
action between the craters. For example, the first crater did 
not collapse as a result of the second impact, nor was the 
first crater erased as a result of shaking or ejecta filling in 
the first crater. 

The authors conclude that large impacts onto compress- 
ible, highly porous targets may not reach the gravity regime 
in which the gravity scaling laws can be employed to predict 
crater diameter and depth. In Fig. [3] the compaction craters 
in perlite and mixtures of sand, perlite and fly ash (open 
symbols) are strength-dominated. As a result, a porous, 
compressible object may have a very high resistance to dis- 
ruption even if both the tensile and compressive strengths 
are low. 

The occurrence of compaction cratering in nature is not 
understood and presently a subject of debate. More work 
on the compaction cratering phenoma is needed. If com- 
paction cratering is prevalent, the bulk density of a porous 
compressible object may be significantly increased over the 
age of the solar system by compaction from impacts. 

5.1.4. Cratering in other volatile materials 

KBOs show wide diversity in volatile content. Large 
KBOs that are bright enough for detailed spectroscopic 
study show evidence of significant volatile content (for ex- 
ample, methane and ethane ices) (Barucci et al. this vol- 
ume, \Brown et a/.||2007a| \Barucci ~ei al. 2005). Labora- 
tory experiments have found that the addition of material 



more volatile than H2O ice, such as CO2 and NH3 (Burchell 



\et al.\l998\\Burchell and Johnson\2005\\Schultz\\996) can 
increase the strength of the target and as a result decrease 
the cratering efficiency. 

The phase of the volatile is also important. Comet nu- 
clei and their precursors may contain trapped pockets of gas 
under high internal pressures. If an impact event releases 
trapped gas, the vapor expansion may aid in the ejection of 
more mass than would be possible from the kinetic energy 
of the impact itself ( Durda et a/.|2003| Schultz et a/.|2005] 



\Holsapple and Housen\2007) . 
5.1.5. Cratering in layered targets 



Belton et al. (2007i suggest that all three Jupiter Fam- 
ily Comet nuclei (believed to originate in the scattered disk 
component of the Kuiper Belt) that have been closely ob- 
served to date (Wild 2, Borrelly, and Temple 1) show evi- 
dence of layering. Belton et al. (2007) propose that this lay- 
ering is primordial and a result of the accretion process. By 
extrapolation, the precursor objects in the Kuiper Belt may 
also have layered structures. Whether the observed layering 
is primordial or not is a matter of debate; however, surface 
layering (a devolatilized "crust") was predicted for comets 
based on thermal evolution models (Belto n~and A'Hearn\ 
1999). Layering of different strength materials does ex- 



plain features seen on other objects in the solar system. For 
example, concentric crater morphology on the moon and 
slightly filled in linear structures on the asteroid Eros can 
be explained by regolith overlying more competent rock. 
Oberbeck and Quaide\§l967\ and\Ryan et al. ( 199ll con- 



ducted experiments on regolith covered targets and deter- 
mined that the morphology of the resulting crater changed 
depending on the depth of the regolith. This result has been 



confirmed with numerical experiments by\Senft and Stewart 



(2006). [fcyaw et a/.| ( |1991| l conducted drop tests to study the 
collision outcome of aggregate projectiles impacting differ- 
ent depths of regolith (fine particles overlaying a concrete 
surface). When the depth of the regolith was at least the size 
of the projectile, the aggregate lost <10% of its mass when 
dropped from a height that resulted in catastophic disrup- 
tion when the surface was not covered with a layer of fine 
particles. The porous regolith was very efficient at dissipat- 
ing the impact energy. 

In addition, impact experiments into granular mixtures 
of H2O ice, CO2 ice, and pyrophylite that have experi- 
enced thermal stratification produce craters with very differ- 
ent morphologies { Arakawa et aZ. |[2000 1. Finally, {Schultz 



2003) looked at the effect of layering on crater scaling. 
Craters retain their original diameter until the layer be- 
comes less than twice the projectile diameter (for vertical 
impacts) or less than a projectile diameter (for oblique im- 
pacts). Even though the final crater depth is limited by 
the substrate, the diameter remains unaffected. Imagery of 
craters on the surfaces of KBOs would provide information 
about near-surface layering. 
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5.2. Catastrophic Disruption 

As mentioned above, there has been much less work 
in the catastrophic disruption regime than the cratering 
regime. We begin this section with a brief summary of 
the laboratory experiments of catastrophic disruption ei- 
ther using ice targets or investigating a range of porosities. 
Next, numerical experiments on ice or porous targets are 
presented. 

5.2.7. Catastrophic Disruption Laboratory Experiments 

Strength-regime laboratory experiments have investi- 
gated the catastrophic disruption of icy and porous targets. 
Both Love et al.\ JT993) and \Ryan et al.\ (119991 conducted 



catastrophic disruption impact experiments into macrop- 
orous targets. Arakawa et al. (2002) performed impact 



experiments into nonporous and porous pure ice and ice 
silicate mixtures. In most of the experiments, the porous 
targets were more difficult to disrupt because the kinetic 
energy of the projectile is partitioned into crushing energy 
to fill void spaces { Love et aT||1993 i and the shock wave 
reflects off of the large number of free surfaces. The result 
is significant attenuation of the shock wave compared to 
solid materials. 



Ryan et al. ( 19991 conducted 20 low speed (100 m s _1 ) 



impact experiments into solid and porous ice targets. They 
found that porous ice targets, though significantly weaker 
than the solid ice targets under static conditions, had a dis- 
ruption strength equivalent to the solid targets with similar 
total mass. The authors attribute this behavior to the effi- 
cient dissipation of energy in void spaces. 



Love et al. ( 1993 i ran a series of hypervelocity experi- 



ments (4.8-6.0 km s _1 ) into glass targets of varying poros- 
ity and strength. They found that the specific energy needed 
to catastrophically disrupt the target was proportional to (1- 
porosity) -3 6 . Their results suggest that the porosity of the 
target is more important for the collision outcome than the 
compressive strength of the target. More work is needed to 
separate the effects of porosity and strength. Impacts into 
the more porous targets resulted in deeper penetration of 
the projectile but the excavated volume was about the same 
as in less porous targets. With higher porosity, the damage 
from the impact was more localized. These results suggest 
that porous objects in the solar system would have longer 
lifetimes against collisional disruption than monoliths of the 
same material. 



periments. However, more work is needed to understand 
how porosity strengthens a material and how to predict dis- 
ruption strength as a function of porosity. 

5.2.2. Catastrophic Disruption Numerical Simulations 

Investigations of catastrophic disruption in the gravity 
regime rely upon numerical experiments. Studies including 
KBO analog materials are limited. Aspha ug et q/.| (fl998 ) 
have considered km s _1 impacts into macroporous targets. 



Benz and Asphaug (1999) and Leinhardt and Stewart, in 



prep, have studied the disruption of solid ice targets. More 
complex simulations of KBOs including microporosity and 
mixed silicates with ice have yet to be conducted. 

Using a SPH code, -Asphaug et al. (1998]) investigated 
how different internal configurations affect the collision 
outcome. They considered 5 km s _1 rocky impacts onto 
a target shaped like asteroid Castalia, which appears to be a 
contact binary. The possible internal structures considered 
were (i) a solid rock, (ii) a global rubble pile with 50% bulk 
porosity, and (iii) two solid rock pieces separated by a zone 
of highly damaged rock. In all three cases, the mass of the 
target was constant (the density of the rock was changed). 
The model included material strength and self-gravity. 

In the rubble pile case, some of the energy generated by 
the impact is partitioned into collapsing void space. In ad- 
dition, the shock wave reflects off of the free surfaces of 
the rubble pieces. As a result, shock effects were focused 
close to the impact site and the shock pressures were dissi- 
pated much more quickly compared to the solid rock target, 



in agreement with laboratory impact experiments 5.2. 1 



The velocities of the ejected material were higher in the rub- 
ble pile configuration than the solid rock case, resulting in 
a small ejecta blanket or none at all. In the two solid piece 
model, the damaged region in the middle of the body re- 
flects the shock wave so that the damage is localized to the 
piece that was impacted. This study elegantly demonstrates 
the importance of internal structure in the outcome of colli- 
sion events. 

A significant problem limiting numerical studies of hy- 
pervelocity catastrophic disruption events is the vast differ- 
ence in dynamical times between the shock propagation and 
gravitational reaccumulation. To make the problem more 
tractable, Leinhardt and Stewart, in prep, have begun using 
a hybidized shock physics - gravity method to study KBO 
analog objects (Fig. [6]). The impact and deformation stage is 



\Arakawa et al.\ (|2002) performed moderate speed im- E^^ed usin S a shock P h y sics code ' CTH ^Glaun et al. 



pact experiments (150-670 m s _1 ) into ice and ice-silicate 
mixtures to quantify the effect of porosity on disruption 
strength. In pure ice targets, the disruption strength in- 
creased with increasing porosity. Puzzlingly, in mixed ma- 
terial targets, the disruption strength decreased with in- 
creasing porosity. These experiments suggest that the na- 
ture of the material bonding (and material strength) can be 
as important as the bulk porosity. 

The work to date demonstrates that porosity plays a sig- 
nificant role in the outcome of catastrophic disruption ex- 



1990 ), and the results are handed to a A-body gravity code, 



pkdgrav \Stadel\2QQ\\ Richards on et aZ.|2 000; Lein hardt\ 
et a/.|200"0) ). This method allows detailed modeling of the 
shock deformation including heating, phase changes, and 
mixing of material as well as the final gravitational reaccu- 
mulation of fragments. 

These simulations record the provenance of the material 
in the largest remnants and track the degree to which the 
reaccumulated material is processed by the initial impact 
event. For example, the peak shock pressure (and hence 
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Fig. 6. — An example of a hybrid hydrocode to iV-body numer- 
ical simulation of the catastrophic disruption of a solid ice target. 
The first row shows the positions of the target (grey) and projec- 
tile (black) at time and 30 s. The second row shows the target 
after handoff from the hydrocode to the iV-body code at 60 s and 
the largest remnant (r = 35 km) at 140 hr. The projectile (r = 8.4 
km) hit the target (r = 50 km) at 45 degrees and 3 km s -1 . The 
first three frames show a slice through the 3D target and projec- 
tile along the y=0 plane. The final frame shows the surface of the 
largest post-collision remnant. The color coding in the iV-body 
frames show the peak pressure attained by each mass element. 
Blue is the lowest peak pressure (2 x 10 6 dynes cm -2 ), and red 
is the highest (1 x 10 11 dynes cm -2 ). The peak pressure is stored 
in Lagrangian tracer particles during the hydrocode component of 
the simulation. The few non-colored (grey or white) particles in 
the iV-body images are blocks from the hydrocode grid that did 
not contain any tracer particles. The surface of the largest remnant 
shows a mixture of high and low levels of shock deformation. 



the amount of melting or vaporization) experienced by each 
mass element is recorded. In a catastrophic impact, a large 
fraction of the original surface is lost, and the original sur- 
face of the target is only maintained at the antipode of the 
largest post-collision remnant. The surface materials on 
the largest remnant reflect heterogeneous shock processing 
(Fig. [6}. Both highly and weakly shocked material lines 
the surface, while the interior material has a more homoge- 
neous history of moderate shock levels. This suggests that 
the surface materials on KBOs that have suffered a catas- 
trophic impacts could be heterogeneously devolatilized in 
comparison to the interior. The surface heterogeneity may 
also lead to color variations. 

6. Summary and Future Directions 

In this chapter, we have discussed the present state of 
knowledge about the possible physical effects of collisions 
in the Kuiper Belt. The body of work on impact cratering 
and disruption events in KBO analog materials demonstrate 
that composition and internal structure (particularly poros- 
ity) have significant affects on the final outcome of collision 
events. Understanding the role of collisions in changing the 
composition and structure of KBOs is important because 
KBOs are the best representatives of the planetesimals that 
accreted into the outer solar system planets. 

The range of possible outcomes of collisions in the 
Kuiper belt region is more complicated than in the asteroid 
belt. In particular, the low relative impact velocities, the low 
mean density of KBOs, and the likely presence of a variety 
of internal structures, are not fully accounted for in present 
impact models. Impact cratering scaling laws and catas- 
trophic disruption criteria that have been developed for hy- 
pervelocity collisions on solid planetary surfaces and within 
the asteroid belt may not be widely applicable to KBOs. For 
example, the unanticipated large amount of mass ejected 
from comet Temple 1 from the Deep Impact mission re- 
vealed that important physical processes are missing from 
the crater scaling laws ^Holsapple and Housen\2007) . 

We close with recommendations for areas of future work 
to advance our knowledge of the properties of Kuiper Belt 
Objects: 

1 . What is the role of porosity in the outcome of crater- 
ing and catastrophic disruption collisions? And how 
do we separate the effects of porosity and strength? 

Predicting the outcome of collisions into porous ma- 
terials of various strengths requires an improved 
understanding of (i) energy coupling into the tar- 
get and (ii) shock-induced damage (degradation of 
strength). As both are difficult to model accurately 
in codes, clever laboratory experiments that include 
direct measurements of shock wave decay, damage, 
and final crater sizes, in targets that vary porosity and 
strength independently are necessary. In cratering 
events, the residual strength of the damaged target 
ultimately determines the final crater shape and the 
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transition between strength and gravity-dominated 
regimes. For disruption events, laboratory results 
will need to be incorporated into numerical simula- 
tions. 

2. What is the correct way to scale laboratory experi- 
ments on porous materials to larger scales? 

Two length scales appear to dominate the problem: 
(i) the length scale of the porosity with respect to 
the deforming shock wave (microporosity vs. macro- 
porosity) and (ii) the depth of energy coupling of 
the projectile. Laboratory and numerical experiments 
can directly address the effects of varying each length 
scale as the problem size increases from laboratory 
targets to planets. More information about the scale 
of porosities in KBOs can be obtained from studies 
of comet nuclei and inferences of rubble pile vs. dif- 
ferentiated internal structures in large KBOs. 

3. How should volatility be incorporated into scaling 
laws? 

Vapor generation (or release) from a collision event 
would affect the momentum of the flow of exca- 
vated or dispersed material. This difficult problem 
needs more information on the actual composition of 
KBOs. 

4. How does differentiation or layering affect the catas- 
trophic disruption threshold? 

The propagation of the impact shock wave through 
the target is influenced by a layered internal structure. 
This tractable problem can be addressed through lab- 
oratory and numerical experiments of plausible inter- 
nal configurations in KBOs. 

5. How does the mass ratio of the colliding bodies 
change the catastrophic disruption criteria in the hy- 
pervelocity regime? 

The shock pressure profile through the target depends 
on the size and velocity of the projectile. Laboratory 
and numerical experiments can directly address this 
problem for solid bodies. Solutions to question 1 in 
this list are required for highly porous bodies. 

6. How can we validate numerical simulations in the 
gravity regime? 

Crater scaling laws have been validated by high grav- 
ity (centrifuge) experiments. In the study of highly 
porous and weak materials, experiments in vacuum 
and under low gravity are also needed. Validation 
of catastrophic disruption simulations in the gravity 
regime will require new techniques. 

7. What is the magnitude of modifications of KBOs 
from mutual collisions compared to other "weather- 
ing" processes? How different are present day KBOs 
from the primordial planetesimals in the outer solar 
system? 



Cumulative changes in observable properties of 
KBOs, including densities, colors, composition, and 
internal structures, can be addressed by updating col- 
lisional evolution models of the Kuiper Belt with 
the latest understanding of collisional processes in 
porous, icy bodies. Given the wide range of possible 
physical properties of KBOs, studies of individual 
collisions are warranted to examine common impact 
scenarios. At present there is no certain answer, and 
our understanding will be driven by observations to 
come. 
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